Search for: All records

Abstract. The Laurentian Great Lakes significantly influence the climate of the Midwest and Northeast United States due to their vast thermal inertia, moisture source potential, and complex heat and moisture flux dynamics. This study presents a newly developed coupled lake–ice–atmosphere (CLIAv1) modeling system for the Great Lakes by coupling the National Aeronautics and Space Administration (NASA) Unified Weather Research and Forecasting (NU-WRF) regional climate model (RCM) with the three-dimensional (3D) Finite Volume Community Ocean Model (FVCOM) and investigates the impact of coupled dynamics on simulations of the Great Lakes' winter climate. By integrating 3D lake hydrodynamics, CLIAv1 demonstrates superior performance in reproducing observed lake surface temperatures (LSTs), ice cover distribution, and the vertical thermal structure of the Great Lakes compared to the NU-WRF model coupled with the default 1D Lake Ice Snow and Sediment Simulator (LISSS). CLIAv1 also enhances the simulation of over-lake atmospheric conditions, including air temperature, wind speed, and sensible and latent heat fluxes, underscoring the importance of resolving complex lake dynamics for reliable regional Earth system projections. More importantly, the key contribution of this study is the identification of critical physical processes that influence lake thermal structure and ice cover – processes that are missed by 1D lake models but effectively resolved by 3D lake models. Through process-oriented numerical experiments, we identify key 3D hydrodynamic processes – ice transport, heat advection, and shear production in turbulence – that explain the superiority of 3D lake models to 1D lake models, particularly in cold season performance and lake–atmosphere interactions. Critically, all three of these processes are dynamically linked to water currents – spatially and temporally evolving flow fields that are structurally absent in 1D models. This study aims to advance our understanding of the physical mechanisms that underlie the fundamental differences between 3D and 1D lake models in simulating key hydrodynamic processes during the winter season, and it offers generalized insights that are not constrained by specific model configurations. 

Extratropical cyclones develop in regions of enhanced baroclinicity and progress along climatological storm tracks. Numerous studies have noted an influence of terrestrial snow cover on atmospheric baroclinicity. However, these studies have less typically examined the role that continental snow cover extent and changes anticipated with anthropogenic climate change have on cyclones’ intensities, trajectories, and precipitation characteristics. Here, we examined how projected future poleward shifts in North American snow extent influence extratropical cyclones. We imposed 10th, 50th, and 90th percentile values of snow retreat between the late 20th and 21st centuries as projected by 14 Coupled Model Intercomparison Project Phase Five (CMIP5) models to alter snow extent underlying 15 historical cold-season cyclones that tracked over the North American Great Plains and were faithfully reproduced in control model cases, providing a comprehensive set of model runs to evaluate hypotheses. Simulations by the Advanced Research version of the Weather Research and Forecast Model (WRF-ARW) were initialized at four days prior to cyclogenesis. Cyclone trajectories moved on average poleward (μ = 27 +/− σ = 17 km) in response to reduced snow extent while the maximum sea-level pressure deepened (μ = −0.48 +/− σ = 0.8 hPa) with greater snow removed. A significant linear correlation was observed between the area of snow removed and mean trajectory deviation (r2 = 0.23), especially in mid-winter (r2 = 0.59), as well as a similar relationship for maximum change in sea-level pressure (r2 = 0.17). Across all simulations, 82% of the perturbed simulation cyclones decreased in average central sea-level pressure (SLP) compared to the corresponding control simulation. Near-surface wind speed increased, as did precipitation, in 86% of cases with a preferred phase change from the solid to liquid state due to warming, although the trends did not correlate with the snow retreat magnitude. Our results, consistent with prior studies noting some role for the enhanced baroclinity of the snow line in modulating storm track and intensity, provide a benchmark to evaluate future snow cover retreat impacts on mid-latitude weather systems.